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ABSTRACT 

In the presence of a strong m = 2 component in a rotating galaxy, the phase space 
structure near corotation is shaped to a large extent by the invariant manifolds of 
the short period family of unstable periodic orbits terminating at Lj or L 2 . The main 
effect of these manifolds is to create robust phase correlations among a number of 
chaotic orbits large enough to support a spiral density wave outside corotation. The 
phenomenon is described theoretically by soliton-like solutions of a Sine-Gordon equa- 
tion. Numerical examples are given in an N-Body simulation of a barred spiral galaxy. 
In these examples, we demonstrate how the projection of unstable manifolds in con- 
figuration space reproduces essentially the entire observed bar-spiral pattern. 
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1 INTRODUCTION 

Any gravitational theory of galactic spiral structure based 
■ on orbits (see Contopoulos 2004, pp. 497-503, for a review) 
must provide a model of the way the phases of the orbits, i.e, 
the angular positions of the apsides, become correlated, so 
' as to reproduce self-consistently the observed spiral pattern. 
Models for normal spirals based on stable periodic orbits of 
the x\ family, with apocenters aligned along the spiral, were 
proposed by Contopoulos and Grosb0l (1986, 1988). In such 
models, the termination of the main spiral is placed near 
the 4/1 resonance (Contopoulos 1985, Patsis et al. 1991, 
1994, Patsis and Kaufmann 1999), while weak extensions 
may also be found beyond the 4/1 resonance, reaching the 
corotation region. On the other hand, such models are prob- 
ably not applicable at all if the spirals are supported mostly 
by chaotic orbits, as, for example, in the case when a bar is 
present (Kaufmann and Contopoulos 1996). The dominance 
of chaotic orbits near and beyond corotation renders these 
orbits 'instrumental' (Kaufmann and Contopoulos 1996) in 
successful self-consistent models of barred-spiral galaxies. In 
particular, chaotic orbits with Jacobi constants exceeding 
the value at Li^ can circulate in and out of corotation, thus 
supporting both the spiral and the bar (e.g. the 'hot' popu- 
lation of stars in the N-Body models of Sparke and Sellwood 
1987). 

A central question in the study of disk galaxies is 
whether the bar drives spiral dynamics (e.g. Goldreich and 
Tremaine 1978, Athanassoula 1980, Schwarz 1984, 1985), or 
the spiral structure is a recurrent collective instability char- 



acterized by its own independent dynamics (coupled per- 
haps to the bar dynamics, e.g. Sellwood 2000 and references 
therein). In order to understand this problem, a relevant 
key question is what mechanisms generate the spiral pattern 
and whether such mechanisms work preferably with regu- 
lar or with chaotic orbits. In the present paper we address 
the question of mechanisms generating phase correlations of 
chaotic orbits able to support a spiral pattern. 

The set of chaotic orbits near corotation fill stochasti- 
cally a connected chaotic domain of the phase space. The 
stochastic character of the orbits could be naively perceived 
as opposed to any persistent phase correlation among these 
orbits. However, it is well known in the theory of dynamical 
systems that the effective randomness exhibited by chaotic 
orbits is to a certain extent only apparent. This is because 
of at least two reasons: 

i) The loss of information along a chaotic orbit becomes 
important only for times longer than the Lyapunov time of 
the orbit (the inverse of the Lyapunov Characteristic Num- 
ber). 

ii) Even for times longer than the Lyapunov times the 
loss of information in a chaotic domain, due to the expo- 
nential stretching of a phase space volume element, which 
is maximum along one particular local direction, is accom- 
panied by a contraction of the element across the same di- 
rection, i.e a gain of information across this direction. Such 
local directions are determined by the invariant manifolds of 
unstable periodic orbits in the chaotic domain. The concept 
of invariant manifolds is a key concept in understanding the 
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chaotic dynamics in the corotation region as will be analyzed 
in section 2 below. 

In particular, we shall provide evidence that the unsta- 
ble invariant manifolds of the short-period unstable periodic 
orbits, which form two families, one terminating at L\ and 
the other at L2 (see section 2), provide a mechanism yielding 
phase correlations of chaotic orbits that support, precisely, 
a spiral density wave beyond the bar. We find that the in- 
tersections of invariant manifolds with configuration space 
define patterns quite similar to the spiral pattern of the ob- 
served density field. By using an analysis parallel to the 
analysis given in Voglis (2003), it is theoretically demon- 
strated (subsection 2.2) that the angular positions of the 
apocenters of chaotic orbits with initial conditions along an 
unstable manifold can be given, as a function of time, by a 
soliton-type solution of a Sine-Gordon equation. This solu- 
tion is quantitatively accurate only for very small amplitudes 
of the perturbation. However, even for large amplitudes of 
the perturbation we find the same qualitative behavior, that 
is the successive apocenters of chaotic orbits remain corre- 
lated for quite long times, despite the positive Lyapunov 
exponents of these orbits. This phase correlation causes the 
manifold to support the spiral pattern. We furthermore pro- 
vide a numerical calculation of invariant manifolds in an 
N-Body simulation of a barred spiral galaxy. In this calcula- 
tion, we superpose the figures of the manifolds on the locus 
of local maxima of the spiral density given directly by the 
N-body particle distribution projected on the plane of the 
disc. These two figures show a satisfactory agreement. 

A general property of invariant manifolds is that the un- 
stable manifold of one periodic orbit cannot intersect itself 
or the unstable manifolds of other periodic orbits. When the 
phase space is compact (i.e. no escapes are possible), the lat- 
ter property implies also that the manifolds of all the orbits 
fill densely the phase space and that they are locally parallel 
to each other (the same property holds also for stable mani- 
folds). We find, with numerical examples, that this geometric 
arrangement of the manifolds is retained to a large extent 
even when the phase space is not compact, i.e., some escapes 
are possible, provided that the rate of escapes is slow. Given 
this property, we anticipate that the effects shown below for 
a particular class of manifolds corresponding to simple peri- 
odic orbits forming families terminating at L\ (or L2) should 
in principle be present also in the manifolds of many other 
families of unstable periodic orbits. Finally, another impor- 
tant property of invariant manifolds is their structural sta- 
bility, i.e., the manifolds retain their phase-space geometric 
structure under small perturbations of a system. 

The paper is organized as follows: In section 2 we de- 
fine and calculate the invariant manifolds of unstable pe- 
riodic orbits in a simple Hamiltonian model describing the 
dynamics in the corotation region. In particular, we discuss 
the role played by invariant manifolds in the dynamics of 
chaotic orbits. Section 3 presents numerical results from an 
N-Body simulation of a barred spiral galaxy. The manifolds 
are calculated first in a Poincare section and then their in- 
tersections with the 2D configuration space of the disc are 
plotted. These plots are compared to the spiral pattern of 
the system as derived by the particle distribution on the 
disc. Section 4 summarizes the main conclusions from the 
present study. 



2 THEORY 

2.1 The phase space structure near corotation 
and the form of invariant manifolds 

The general form of the Hamiltonian of motion in the rotat- 
ing frame of a disc galaxy with a non-axisymmetric potential 
perturbation is: 

H(r,e,p r , Pe ) = i(p, 2 + jf ) - n pP e + V (r) + Vi(r,0) (1) 

where (r,9,p r ,pg) are polar coordinates and their respective 
conjugate canonical momenta, and Q p is the angular velocity 
of the pattern. Assuming a bar-like perturbation V\ , a good 
first approximation to the Hamiltonian ([lj in the neighbor- 
hood of corotation is given by (Contopoulos 1978): 

H = k*Ii + a, 1 1 + 2b*hh + da + A, cos 2<9 2 (2) 

In this formula, 7i is the action corresponding to epicyclic 
oscillations, ft* being the respective epicyclic frequency. I2 
is an 0(Aa.) correction to the quantity pg — J*, where pg is 
the angular momentum of an orbit in the inertial frame and 
J, is the value of pg for a circular orbit at the corotation 
radius of the unperturbed (Vi = 0) potential. The angle 02 
is also an 0(^4*) correction to the azimuthal angle 0. For 
simplicity, in the analysis below we set I2 ~ pe — J* and 

0a » e. 

Since the epicyclic angle 0i, conjugate to I\, is ignorable 
in the Hamiltonian ((2]), this Hamiltonian model is integrable, 
7i being a second integral independent of and in involution 
with 77. However, if any other term exp(i(n0i + 17162)) of 
the harmonic expansion of Vi is added to the Hamiltonian 
(O, the Hamiltonian becomes, in general, non-integrable and 
chaos is introduced to the model. As an example, we may 
consider the Fourier modes n = l,m = ±2, which are im- 
portant near Lindblad resonances. A simple example of a 
Hamiltonian model including such terms is 

H = k*Ii + a,ll + 2b,hl 2 + e, l\ + A, cos 2<9 2 

-e( 2(Jl + /lo) ) X/2 l2[cos(6>i + 2&) - cos(0! - 20 2 )] (3) 

The dependence of the Fourier terms cos(0i ±282) on the ac- 
tions 7i , I2 introduced in Eq. Q is a simplified model used to 
demonstrate a number of phenomena relevant to the discus- 
sion below, while, in reality, the dependence of any Fourier 
term of Vi on the actions would be determined by the partic- 
ular functional form of Vi . iio is just a smoothing constant 
to avoid the singularity of the derivative with respect to 7i. 
At any rate, the crucial parameter in all these phenomena 
is the overall size of the perturbation terms cos(0i ± 202), 
determined in Eq.@ by e. 

For Jacobi constant values larger than the value at the 
unstable Lagrangian equilibrium points Li,a, there is a fam- 
ily of unstable short-period orbits around L\, or L2 (Fig. 1, 
schematic) . These orbits form loops of size that increases as 
the value of the Jacobi constant increases. We denote with 
7^1,2 the point where a short-period orbit passes from its 
apocenter. We shall also refer to the orbits themselves as 
the short period orbits Pl\ (or Pm). The curves with ar- 
rows correspond to the equipotential passing from 71,2 and 
the arrows indicate the sense of flow of orbits inside corota- 
tion (prograde, counterclockwise) or outside corotation (ret- 
rograde, clockwise). A similar family of stable short-period 
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orbits with apocenters Pl4,5 exists also for the stable La- 
grangian points Z/4,5 for values of the Jacobi constant ex- 
ceeding the value at L4.5. 

Now, a visualization of the phase space of the system 
Q can be obtained by means of a suitable Poincare surface 
of section. Figure 2a shows the surface of section (12,62) 
corresponding to the successive crossings of the angle 81 
along an orbit with the values 2n7r, n = 1,2,..., i.e., at 
the apocenters of orbits, when e is very small (e = 10 -3 ). 
The basic phase portrait (Fig. 2a) resembles a pendulum 
phase portrait, with a separatrix dividing the phase space 
into libration and rotation regions. The fixed points Pli, 2 
and Pl4,5 correspond to the same points as in Fig. 1. The 
points Pli,2 exist for a range of values of the Jacobi con- 
stant Ej > Ej.li,2, while the points Pl4,5 exist for a range 
of values of Ej, such that Ej > Ej,l4,,5 > Ej,li,2- In the 
latter case, the phase portrait is divided in two libration 
domains, around the stable points Pl4,5, and two rotation 
domains, above and below the unstable fixed points Pli,2, 
in which the orbits make rotations around the galaxy in ei- 
ther the retrograde sense (above Pli, outside corotation), or 
prograde sense (below Pli, inside corotation). On the other 
hand, if Pj,l4,5 > Ej > Ej,li,2, there is an energetically 
forbidden domain that replaces some of the KAM curves 
around Pla and Pls of Fig. 2a. 

Figure 2a gives initially the impression of being com- 
posed only by invariant curves, i.e., by regular orbits. How- 
ever, a closer look reveals that there is a small degree of chaos 
in the neighborhood of the unstable points Pli, 2- This chaos 
is clearly manifested by computing the form of the invariant 
manifolds emanating from Pli or Pl2 (Fig. 2b). 

By definition, we call unstable manifold W u of the short 
period unstable periodic orbit Pli the set of points in phase- 
space with pre-images approaching asymptotically the orbit 
Pli as t — > —00. Similarly, the stable manifold W s of Pli 
is the set of points in phase-space with images approaching 
asymptotically the orbit Pli as i -> 00. Both sets, W" 
and W s are two-dimensional manifolds embedded in the 
three-dimensional Jacobi constant hypersurface of the four- 
dimensional phase-space. 

If, now, we work on a particular Poincare surface of sec- 
tion, such as, for example, the section of the apocenters as 
in fig. 2, the unstable periodic orbit Pli appears as an unsta- 
ble fixed point Pli - Furthermore, the two-dimensional mani- 
folds W u or W s intersect this section along one-dimensional 
submanifolds on the two-dimensional surface of section. 
These submanifolds are the set of points approaching asymp- 
totically the fixed point Pli as t — > —00 or t — > 00, respec- 
tively. In order to distinguish between two-dimensional man- 
ifolds in phase space and one-dimensional manifolds in the 
surface of section, we call the latter apocentric or pericen- 
tric manifolds, depending on whether we take the surface of 
section at the apocenters or pericenters of the orbits respec- 
tively. In fact, one can define a surface of section in many 
different ways, but the definitions adopted here, in terms of 
the apocenters or the pericenters, are more convenient for 
the present study. 

A well known theorem of dynamics (Grobmann 1959, 
Hartman 1960) proves that: 

a) in autonomous Hamiltonian systems the manifolds 
on the Poincare surface of section are invariant, i.e., coincide 
with their images under the Poincare map, and 



b) they approach the fixed point Pli in the directions 
tangent to the eigenvectors of the linearized map around 
Pli- In particular, the unstable manifold is tangent to the 
eigenvectors associated with the absolutely larger real eigen- 
value (say Ai) of the monodromy matrix at Pli, while the 
stable manifold is tangent to the eigenvectors associated 
with the absolutely smaller eigenvalue, A2, of the same ma- 
trix. By the symplecticity property of Poincare maps we 
furthermore have A1A2 = 1, thus |Ai| > 1 and | A2 1 < 1. 

Property (b) can be exploited to calculate numerically 
the invariant manifolds emanating from a particular unsta- 
ble periodic orbit. Namely, let the functions 

8' 2 = F(8 2 ,h), I' 2 = G(6 2 ,l2) (4) 

represent the Poincare map of the system under study (the 
functions F, G, are determined numerically). We then cal- 
culate, by a Newton-Raphson scheme, the position of one 
unstable periodic orbit, say Pli, which corresponds to one 
root (#2,0,^2,0) °f the set of equations 82,0 = F(8 2 fi, I 2 ,o), 
^2,0 = G(82 t o,l2,o)- The monodromy matrix at the fixed 
point Pli is defined as 

/ dF_ dF_ \ 

A= H M (5) 

\ d0 2 dl 2 / (0 2 O ,I 2 O ) 

where partial derivatives can also be determined numeri- 
cally. The eigenvalues Ai and A2, of A, are real and recip- 
rocal. We then take a small initial segment dS with one of 
its endpoints coinciding with the fixed point Pli and the 
other endpoint oriented on the direction of the eigenvector 
associated with one of the eigenvalues Ai, or A2. Finally, we 
consider a large number of initial conditions of orbits along 
the segment dS. The successive Poincare map images of the 
points along the segment associated with Ai numerically de- 
termine the unstable one-dimensional invariant manifold of 
Pli- Similarly, pre-images of points along the segment asso- 
ciated with A2 determine the stable one-dimensional invari- 
ant manifold of Pli. 

Figure 2b shows a focus near Pli, 2 of the apocentric in- 
variant manifolds emanating from Pli, 2 for the same model 
and parameter values as in Fig. 2a. The form of the mani- 
folds of Pl2, when they approach the neighborhood of Pli, 
is characteristic of what is known as 'homoclinic chaos' in 
the theory of Hamiltonian dynamical systems. Namely, the 
stable and unstable manifolds have transverse intersections, 
called homoclinic points. The part of a manifold between any 
two successive homoclinic points determines a lobe. The suc- 
cessive Poincare images of the lobes are more and more elon- 
gated along the two eigendirections, creating a very compli- 
cated homoclinic tangle. That the complexity of the dynam- 
ics induced by the homoclinic tangle is the source of chaos in 
Hamiltonian systems was already pointed out by Poincare 
(1892). However, a non-schematic numerical calculation of 
invariant manifolds in a simple Hamiltonian model was given 
relatively recently (Contopoulos and Polymilis 1993). It has 
been established theoretically that this dynamics is isomor- 
phic to the dynamics of the so-called A— set (e.g. Contopou- 
los 2004, p. 152), which obeys the Bernoulli shift map. The 
latter is considered as a paradigm of chaotic system. 

The geometric structure of invariant manifolds in phase 
space is preserved even for large values of the nonlinearity 
parameter. This effect is shown in Figs. 2c,d. The parameter 
e has now a value 10 -2 , i.e., ten times larger than in Figs. 
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2a,b. As a result, the phase space near Pli is now largely 
chaotic, and most rotational invariant curves are destroyed 
(Fig. 2c). However, the apocentric invariant manifolds em- 
anating from Pli, 2 (Fig- 2d) preserve a geometric structure 
similar to that of Fig. 2b, except for the fact that the overall 
vertical scale of Fig. 2d is larger by a factor five from that 
of Fig. 2b. This implies that, as e increases, the size of the 
lobes of the invariant manifolds increases. The lobes extend 
to bigger distances away from Pli. 2. 

Now, the unstable invariant manifolds of an unstable 
periodic orbit, starting at the directions (U,UU), establish a 
preferential direction of chaotic stretching in phase space in 
the forward sense of time. In particular, all the orbits with 
initial conditions inside the area defined by the lobes of one 
unstable manifold are forced to follow the same preferential 
direction as that of the manifold. Furthermore, the unsta- 
ble manifolds of other unstable periodic orbits in the same 
chaotic domain cannot intersect the unstable manifold of the 
first orbit, thus they also follow the same directions. Such 
directions are defined in the neighborhood of every point in 
the chaotic domain and they serve as the basis of construc- 
tion of numerical indicators distinguishing between ordered 
and chaotic orbits (Voglis et al. 1998, 1999, 2002). 

Returning to the discussion of manifolds emanating 
from Pli, the relevant remark regarding the chaotic behav- 
ior of orbits near corotation is that, as the lobes of the unsta- 
ble manifold U, or UU, approach close to Pli, they become 
more and more elongated, so that they repeatedly cross the 
straight lines corresponding to the directions S, or SS (Figs. 
2b,d). At every such crossing, a lobe goes from the chaotic 
domain above Pl2 (outside corotation), to the domain below 
Pl2 (inside corotation), or vice versa. As a result, all chaotic 
orbits with initial conditions on, or inside, a lobe repeatedly 
go inside and outside of the corotation region. 

The sequence of these events has a Bernoulli shift dy- 
namics, i.e., it is a very chaotic phenomenon. This fact 
notwithstanding, we show in the next subsection that the 
same dynamics imposed by unstable invariant manifolds im- 
plies that the phases of chaotic orbits near an unstable man- 
ifold are organized along the dominant directions of stretch- 
ing, due to the area preservation, and they become well cor- 
related. This effect, precisely, serves as the basis of support 
of a spiral density wave composed by chaotic orbits. 

2.2 Soliton-like motions and phase correlations of 
chaotic orbits 

Voglis (2003) has shown that in the neighborhood of any 
resonance in a disc galaxy, the phases of stars with initial 
conditions close to the separatrix associated with the reso- 
nance follow a soliton-type motion. A particular numerical 
application was given in the case of the Inner Lindblad Res- 
onance, but exactly the same formalism can be applied to 
other resonances. In this section we shall assume the Hamil- 
tonian to be of the form ([3]) and describe soliton-type cor- 
relations applicable to the orbits near Pli, 2- 

To this end, we introduce the canonical transformation 
(0i,02,/i,/ 2 ) -» (0i,0£,Jj,i2) defined by 

' - 2Z 

1 BE ' 



_ dF 



1,2 



(6) 



way as to eliminate the dependence of the Hamiltonian on 
the angle 0i up to terms of order O(e). Precisely, we have 

F = 9 1 l' 1 + e 2 I 2 + e( 2{Ii+Ilo) )^x (7) 

[sin(0i + 20 2 ) - sin(0i - 20 2 )] 

Substituting the transformation ([6| into the Hamiltonian 
the Hamiltonian takes the form: 

H = K#l[ + a*Ii + 2b,l' 1 l' 2 + c/2 2 + A, cos 20 2 + 0(Ae)(8) 

where primed variables differ by 0(e) terms from their re- 
spective non-primed variables, and both quantities A and e 
are assumed small. Ignoring terms of order O(Ae), Hamil- 
ton's equations of motion in the Hamiltonian ((H) can be cast 
in the form: 



ip — luq sin ip = 



(9) 



where ip = 2(0 2 — it/2) and u)q = Sc^A* (the constant c» is 
negative, while A* is positive) . The first integral of (0 reads 



— + u cos ip 



c 



(10) 



The quantity (|10p represents an approximate integral of mo- 
tion which, in a strict sense, is valid only for orbits lying on 
invariant curves such as those of Fig. 2a. In particular, the 
value of the constant C = ojq corresponds to a theoretical 
separatrix curve passing from Pli, 2- Along this curve, the 
solution of ((9(1 reads 

tan (p = exp[±(a;ot - fco£)] (11) 

where the parameter £ is related to the initial phase ipo of 
an orbit on the separatrix via £ = — lntan(?/>o/4), and 
ko = Wo/wi, with uji = k* + 0(A*). This choice of uii allows 
us to follow the evolution of the phases ip for the successive 
Poincare consequents of the orbits when 6\ — uj\t = 2nn, 
with n = 1, 2, . . ., i.e., at the apocenters of the orbits. 

Now, the solution (111[) can be viewed as a kink, or anti- 
kink, soliton-type solution of the Sine-Gordon equation: 



d 2 i> 
dt 2 1 d$ 2 



2& 2 ip _ 2 



2lo sin ip = 



(12) 



where the generating function F(8i,I'i) is chosen in such a 



In Eq. (|12[) , £ was attributed the status of an independent 
variable of the problem expressing the initial value of 02 
on the Poincare surface of section. The reason for this is 
that the number of stars in a galaxy is large enough so as 
to expect that the separatrix defined by the apocenters of 
stars is populated by stars all along its length. We may thus 
view Eq. (|12[) as describing the evolution of the phases ip of 
a chain of stars populating the separatrix. This is given as a 
function of the time, t, but also of the initial phase of each 
star, given by £. 

The above analysis was based on a canonical transfor- 
mation of first order in the perturation e, which also ap- 
proximates parts of the manifolds as separatrices. The ques- 
tion then arises whether the soliton-type flow of phases in- 
duced by Eq. (|12|1 is a good approximation of the true flow 
of phases along the invariant manifolds as calculated nu- 
merically, i.e., without any approximation. In order to ex- 
plore this question, we first avoid the multiplicity of roots 
of Eg ■ P with respect to ip by differentiating this equation 
with respect to time, that is ip = ±2u>osech(u)ot — fco£), or 
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8' 2 — ±u)Qsech{uiot — ko() (13) 

In order to compare numerical results with analytical predic- 
tions, we then plot the theoretical relation (|13p against the 
numerical evaluation of 82 (by Hamilton's equations from 
the Hamiltonian ©) for an ensemble of points being ini- 
tially in a small segment along the unstable eigenvector 
pointing in the UU direction, from Pli, and the U direc- 
tion, from Pl2- As shown in Fig. 3, for a relatively large 
value of the perturbation parameter e = 10 -2 , and at differ- 
ent time snapshots, the numerical solution 8 (displayed by 
dots in Fig. 3), as function of (£,£) follows approximately 
the analytical solution (Eq. (|13[l . displayed by solid line in 
Fig. 3). In fact, the theoretical solution for 62 (Eq[TTJ with 
ip — 29' 2 — tv ~ 202 — vr has four non-communicating branches 
which correspond to the non-communicating parts of a the- 
oretical separatrix that would pass through Pli, 2 (Fig. 2a) 
if chaos was absent, namely a) upper and b) lower part for 
< 82 < 7r/2, and 37r/2 < 62 < 2ir and c) upper and d) 
lower branch for n/2 < 62 < 3tt/2. These are merged to 
two branches for the solution of 82 via Eg . (|13fl . On the con- 
trary, the invariant manifolds allow communication of the 
branches via the expanding lobes as those shown in Fig. 2 
b, d. This means that when the orbits approach very close 
to Pli, 2 they are deflected to a different branch instead of 
approaching asymptotically, as t — > 00, Pli or Pl2- This 
means that at, any time t, the motions of the manifold points 
can be followed only within some limits £m»n ^ £ ^ (.max 
with (max — (min large, but not infinite as in the integrable 
case. Since particles jump to a different branch of the equa- 
tion (13), than the branch they were initially, their phases 
82 are mixed in time. In order to circumvent this problem, 
we redefine the value of ( for each orbit whenever the or- 
bit changes branch, according to Eq.([TTJ, with tp = 282 — 7r 
and 82 equal to its value at the particular moment when the 
particle changes branch. 

The solitary motions shown in Fig. 3 correspond to a 
correlated motion of the successive apocenters of chaotic or- 
bits with initial conditions along the manifold. This phe- 
nomenon, along with the physical meaning of the changes of 
branches, will be discussed in subsection 2.3 below. Further- 
more, we find numerically that this correlation persists to 
a large extent even for relatively large values of the pertur- 
bation e as a consequence of the structural stability of the 
invariant manifolds. It is worth examining the consequences 
of persisting phase correlations for the dynamics induced by 
chaotic orbits in the corotation region. 



2.3 Mapping of the invariant manifolds in 

configuration space and the form of spiral 
arms 

As described in subsection 2.1, the apocentric or pericen- 
tric manifolds are one-dimensional objects embedded in a 
three dimensional space of constant energy. This implies that 
they have one-dimensional projections in any plane of this 
space. In the sequel, besides the surface of section (82,12), or 
(8,pe), we shall consider the one-dimensional projections of 
manifolds in the physical, or configuration space represented, 
e.g., by the canonical position variables (r,8), or x = rcosd, 
y — rsind. This allows to compare the topological charac- 



teristics of the manifolds to morphological characteristics of 
the simulated galaxy. 

In order to understand the geometric connection of the 
surface of section to the configuration space, a key remark is 
that, in the corotation region, the surface of section (I2, 82) 
and the configuration plane (r,8), or x — rcosd, y — rsind, 
are nearly isomorphic. This isomorphy is shown schemati- 
cally in Fig. 4. The domains denoted A,B,C,D in Fig. 4a 
are mapped to the domains A',B',C', and D', respectively, 
in Fig. 4b. The mapping is derived as follows: If the origin 
of epicyclic angles 8\ is selected in such a way so that the 
Poincare section points (12,82) for 81 — 2nir, n — 1,2,... 
correspond to successive apocentric passages of an orbit, we 
have r = at every section point. Replacing pg ~ I2 + J*, 
8 ~ 82 in the general form of the Hamiltonian (Eq.([l|), the 
equation of the Jacobi constant H(r,8,p r — 0,po) = Ej 
may be solved for the variable r, yielding r at the apocen- 
ters as a function of 8,pg (or I2, 82) and of the constant Ej. 
One readily finds that for potential functions correspond- 
ing to the usual form of the rotation curve expected in a 
barred galaxy, the apocentric radius r is a monotonically in- 
creasing function of pg both inside and outside corotation. 
In the inner parts of the galaxy we expect that the velocity 
of circular orbits raises as v c oc r thus r cx p\ . On the 
other hand, beyond corotation the rotation curve is nearly 
flat up to a distant radius, i.e., v c ~ const, yielding r cc pg. 
Finally, in the Keplerian limit v c oc r -1 ^ 2 , thus r oc p\. 
These relations are not seriously distorted by the inclusion 
of epicyclic oscillations. We shall see in the next section that 
these monotonic relations are verified by exact calculations 
in an N-Body system. 

Because of the monotonic relation of r and pg, phase 
portraits in the space (pg,8), or (12,82), can be mapped to 
nearly-isomorphic portraits in the space (r,8), where r, cal- 
culated as indicated above, always represents a local apoc- 
entric distance of an orbit. 

Figures 4a, b show the correspondence between the di- 
rections of the stable and unstable manifolds emanating 
from Pli, 2 in the section (12,82) (Fig. 4a) and the same 
directions in configuration space (x,y) (Fig. 4b). In partic- 
ular, the unstable manifolds U, pointing upwards in Fig. 
4a, are mapped to the directions U of Fig. 4b. In the sep- 
aratrix limit, the unstable manifold U joins smoothly the 
stable manifold S, while the unstable manifold UU joins 
smoothly the stable manifold SS. Furthermore, the branch 
US does not communicate with the branch UUSS. We em- 
phasize that such a theoretical separatrix can be constructed 
for any value of the Jacobi constant at which the periodic 
orbit Pli exists, i.e., the separatrix of Fig. 4b should not be 
confused with an equipotential passing through the equilib- 
rium point L\. The separatrix gives the slow drift of apoc- 
enters along the branch US, or UUSS, that implies a slow 
drift of the guiding centers of epicyclic motion for all orbits 
with initial conditions on these two branches. These separa- 
trices do not coincide with the equipotential curves passing 
through Li, L2 even in the limiting case when Pli, Pl2 tend 
to Li, Z/2. In this respect, the manifolds shown in Fig. 4b are 
in a sense parallel to the 'critical ergos curves' of Lynden- 
Bell and Barot (2003), but at a higher value of the Jacobi 
integral. By an elementary analysis of the linearized flow 
in the neighborhood of the Pli, 2 orbits (named 'Lyapunov 
orbits'), Romero-Gomez et al. (2006) gave plots similar to 
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Fig. 4b. In view of the geometry of the separatrix, these 
authors identify the manifolds as responsible for creating a 
ring structure in barred galaxies like NGC1326. They also 
point out the robustness of these objects, called by them 
'flux rings'. 

In fact, as shown in the previous subsection, the robust- 
ness of the phase flow along a manifold is a consequence of 
the soliton-type motion of the phases of a chain of orbits 
populating the manifold. This is shown in Figs. 4c, d, and 
4e,f, which show the invariant manifolds of Fig. 2d in phase 
space and in configuration space after 4 iterations of an ini- 
tial segment dS = 10~ 4 , (Figs. 4c, d) and after 12 iterations 
(Figs. 4e,f). The correspondence in the form of manifold 
lobes in configuration and phase space is clearly visible in 
these figures. 

We emphasize that each point of figures 4d,f corre- 
sponds to a position of apocenter of a chaotic orbit with ini- 
tial conditions along the manifolds. These positions change 
in time, for each particular orbit, following the soliton-type 
flow described in subsection 2.2. However, the pattern de- 
fined by the manifold remains time-invariant. Furthermore, 
as the lobe oscillations become of larger and larger ampli- 
tude, along the successive iterations of the initial dS, the 
manifold moves far away from the separatrix. Given that the 
sense of rotation of the galaxy is counter-clockwise, we see 
that the outward branches U of the unstable manifolds em- 
anating from Pli,2 define trailing spiral arms beyond Pl\,i- 
This phenomenon will be studied in detail in the next sec- 
tion. In the same time, the inward directions UU support 
the bar, i.e., the same chaotic orbits with initial conditions 
along the manifolds support both the spiral pattern and the 
bar. 

The opening of the manifold lobes is larger and faster 
when e increases. We shall see however, in the next section, 
that the near-isomorphy between phase and configuration 
portraits is preserved for quite large perturbations. In par- 
ticular, in order to study how far can extend the manifold 
spirals, we present a numerical example referring to an N- 
Body simulation of a barred spiral galaxy. It is then possi- 
ble to show that the invariant manifolds U generate patterns 
that closely follow the form of the N-Body spiral pattern of 
the system. 



3 NUMERICAL SIMULATIONS: 
COMPARISON OF INVARIANT 
MANIFOLDS AND SPIRAL ARMS 

Figure 5a shows one characteristic snapshot (for t — 47 in 
units of the half-mass crossing time) of an N-Body simu- 
lation of a rotating barred - spiral galaxy. A series of such 
simulations with N = 130000, is presented in Voglis et al. 
(2006). These experiments are produced with a smooth po- 
tential field code which is an improved version of the code 
of Allen et al. (1990). The initial conditions for a rotat- 
ing system are produced from the end-state of the so-called 
Q— model in the simulations of Voglis et al. (2002). This is 
initially a triaxial configuration produced by a collapse pro- 
cess with cosmological initial conditions consistent with the 
standard A— CDM model. The initial model has an almost 
zero spin parameter A (Peebles 1969). In order to increase 
the total angular momentum of the model, while respecting 



its energy and the scalar virial theorem, the following veloc- 
ity re-orientation process is implemented: At a given snap- 
shot of N-body evolution of the Q- model (tg = 100), the ve- 
locity component v yz of each particle, where the plane (y, z) 
is the plane of intermediate - longest axes, is re-oriented so 
as to become perpendicular to the current cylindrical ra- 
dius r yz of the same particle. After the re-orientation, the 
velocity components v yz of all particles point to the direc- 
tion of rotation clockwise with respect to the x— axis, which 
thus becomes the rotational axis of the system. This is the 
way plots are presented in Voglis et al. 2006. In the plots 
below, however, like Fig. 5a we have reversed the sense of 
the x-axis so that the pattern in these figures is shown to 
rotate counterclockwise (this is the most frequent choice in 
the literature when discussing bar or spiral dynamics). 

The so-produced system, called QR1, is let to evolve 
by the N-Body code until a time tQm = 20 (in units of 
the half- mass crossing time) , and then the same velocity re- 
orientation process is implemented. The new system, QR2, 
is left to evolve again by the N-Body code until tQR2 = 20, 
when the velocity re-orientation process is implemented once 
time. This process can yield different rotating systems which 
have the same binding energy, and nearly the same virial 
equilibrium, but they have consecutively higher and higher 
amounts of angular momentum. 

A detailed analysis of the regular and chaotic orbital 
content as well as of the evolution of the above systems is 
presented elsewhere (Voglis et al. 2006). We only mention 
here the results relevant to the analysis below, which refer 
to the long-time evolution of the QR3 experiment, i.e., the 
experiment with initial conditions provided by three consec- 
utive applications of the velocity re-orientation algorithm: 

a) The initially triaxial configuration has developed a 
thick disc structure, with a thickness ratio of the order of 
0.2. The projection of the system on the disc plane (e.g. Fig. 
5a) clearly shows a bar as well as a grand design spiral pat- 
tern. In most snapshots the calculated spiral pattern speed 
has nearly the same value as the bar speed. The spiral am- 
plitude, however, is variable, ranging from a maximum 2 or 
3 times smaller than the amplitude of the bar to a minimum 
that is nevertheless always above the threshold of statistical 
noise. 

b) In some snapshots the spiral looks detached from the 
bar. Most snapshots, however, manifest a continuity of the 
pattern between the bar and the spiral. 

c) The orbits of all the particles are analyzed as regards 
their regular, or chaotic character, by a combination of effi- 
cient numerical chaos indicators (Voglis et al. 2002). Let us 
note that the orbital analysis of a particular snapshot is pos- 
sible because the code provides an expansion of the potential 
in terms of a set of basis functions, making thereby possible 
the explicit calculation of orbits, variational equations etc. 
The main results are: i) The mass in chaotic motion is found 
to be on the level of 60% to 65% of the total mass, and ii) 
when particles in chaotic or ordered orbits are plotted sep- 
arately, the spiral arms are found to be composed practically 
only by chaotic orbits. On the contrary, the bar is composed 
by both regular and chaotic orbits. 

At the snapshot t — 47, the bar-spiral structure is 
clearly seen by simple visual inspection (Fig. 5a). All dis- 
tances in this figure were divided by the half mass radius 
of the system, which is taken as the unit of length. The 
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thick dots show the positions of local maxima of the pro- 
jected density on the (y, z) plane. This is done by splitting 
the distribution of points in concentric rings r, r + Ar, where 
Ar = 0.2, and calculating the angles of local density maxima 
within each ring. 

Figure 5b shows the distribution of the values of the Ja- 
cobi constant for the set of particles which are located within 
small projected area elements rArAd, with Ad — 2ix/26 
around the values (r, 6) corresponding to the thick dots of 
Fig. 5a. The energy units are derived by the numerical code 
length and time units, i.e., setting the mass of each parti- 
cle equal to m p = 1. In these units, the half-mass cross- 
ing time is thmct = 1.1085 x 10~ 4 , the half-mass radius is 
R h = 9.26 x 10~ 2 . The root -mean square velocity is found to 
be v rm3 = 732.6. The pattern speed (estimated numerically 
by the rate of angular displacement of the density maxima 
of the pattern), is equal to fl p = 27r/9.43 in units of l/t hmct , 
i.e., the period of pattern rotation is about ten times larger 
than the half-mass crossing time (which corresponds to a 
typical period of the particles' orbits). This means that the 
rotating pattern can be safely characterized as an evolving 
density wave. 

Another approximation in the calculation of Fig. 5b is 
the fact that the kinetic energy of a particle includes a v\ /2 
term, which corresponds to the kinetic energy of the velocity 
component perpendicularly to the disc. This small term (a 
few percent of the total kinetic energy) is substracted from 
the evaluation of the Jacobi constant, since, for most par- 
ticles, the motion perpendicularly to the disc is effectively 
decoupled from the motion on the disc plane. We thus ap- 
proximate the Jacobi constant for each particle by: 

E J = \( v l + v 'z) - top(zvy - yv z ) + V(x,y,z) (14) 

where the function V(x,y, z) is calculated by the series of 
basis functions of the N-Body code. 

Figure 5b shows the distribution of Jacobi constants, 
with the above conventions, for the particles located in the 
maxima shown in Fig. 5a, where, in addition, we only count 
the particles in rings beyond the value r = 1, marking, ap- 
proximately, the end of the bar. Thus, the distribution of 
Fig. 5b represents essentially the distribution of values of the 
Jacobi constant for particles located along the spiral arms. 
The peak of this distribution is at a value Ej — — 1.12 x 10 6 , 
which is slightly above the value, -Ej,z,i = —1.133 x 10 6 . 

As described in subsection 2.1, the apocentric manifolds 
are one- dimensional manifolds in the three-dimensional sub- 
space (r,6,pe) of the entire four- dimensional phase space of 
the system. Figure 6 shows a calculation of the unstable 
apocentric manifolds emanating from the periodic orbit Pl2 
in the potential of the system of Fig. 5a (with initial con- 
ditions 6 = 3.2247, pe = 1.194). Three different projections 
of these manifolds are shown, namely a) on the surface of 
section (8,pe), r — 0, r < (Fig. 6a), b) on the plane (r,p$) 
(Fig. 6b) and c) on the configuration space (r,0) (Fig. 6c). 
The Jacobi constant has the value Ej — —1.125 x 10 6 , close 
to the maximum of the distribution of Fig. 5b. (We only 
examine theoretical orbits with x = 0, i.e., lying on the disc 
plane). Fig. 6 shows the seventh Poincare map image of two 
initial segments of length dS = 10~ 4 on the Poincare section, 
tangent to the unstable directions U and UU respectively, 
at the fixed point Pl2- As a property of the manifold, the 



seventh image of the initial dS contains also all the previ- 
ous images from the first to the seventh. This choice of ini- 
tial length was made by the criterion that the numerically 
calculated direction of the segment after its first iteration 
coincides with the eigendirection as calculated by the mon- 
odromy matrix of the periodic orbit Pl2 up to 3 significant 
digits. Furthermore, each segment is populated with 10000 
equally spaced points of initial conditions. Figure 6c shows 
the calculated projections of the apocentric manifolds U and 
UU (2x10000 points) as well as their symmetric manifolds 
with respect to the center y = z = 0, which are the mani- 
folds emanating from Pli- 

Figure 6a is different from the simple picture of mani- 
folds shown in Figs. 2b, d in many aspects. A first difference 
is that the 6 value of Pli is close to the value 8 = it, i.e., 
shifted by 7r/2 with respect to Figs. 2b,d. This is because at 
the snapshot t = 47, the bar is oriented almost horizontally 
(Fig. 5a). But the most serious difference is that the mani- 
folds of Fig. 6a are broken into disjoint pieces. This breaking 
of the image of the manifold on the surface of section is an 
artifact caused by the particular choice of surface of section. 
This can be exemplified with the help of Figs. 7a,b,c,d which 
refer to the same manifold, but calculated for six iterations 
of the initial dS, in order to obtain a clearer figure. Figures 
7c, d show the time evolution of two nearby orbits one of 
which gives the last point, P, of a disjoint piece of the man- 
ifold of Fig. 7a, while the other orbit gives the first point P' 
of the next disjoint piece of the manifold in the same figure. 
We notice that while the orbits are really nearby (Fig. 7d), 
the first orbit has a local maximum of r(t) at the point P, 
while the second orbit has no local maximum of r(t) at times 
near the time of the P— maximum. In fact, the second orbit 
yields a local maximum much later, at the point P' (Fig. 
7d), thus breaking the continuity of Poincare intersections 
in Fig. 7a between the points P and P'. The same phe- 
nomenon is also shown in Fig. 7c, where we clearly see that, 
while the first orbit on the surface of section (r = 0,r < 0), 
represented by the line r = in Fig. 7c, intersects the line 
r = at P, the second orbit passes close to the section line, 
but without intersecting this line. 

We emphasize that the above effect is an artifact of the 
particular choice of surface of section, which is not a true 
Poincare section (Poincare 1892). In fact, it is hardly possi- 
ble to find any surface of phase space in the above model that 
can be proven to be a true Poincare section. Even so, the 
invariant manifold surfaces are continuous in phase space, 
and it is only their intersections with the particular surface 
of section which are discontinuous. Furthermore, the basic 
topological property of invariant manifolds, namely the fact 
that they cannot intersect themselves, is also preserved in 
the surface of section images of the manifold shown in Figs. 
6a or 7a. This fact implies that, as the number of iterations 
increases, i.e., when the manifold is calculated for longer 
times, the manifolds fill the connected chaotic domain em- 
bedding the fixed points Pli, 2 in the surface of section. This 
tendency is clearly observed by comparing Fig. 6a, which de- 
picts the manifold after only one more iteration than Fig. 7a. 
In Fig. 6a the manifold crosses a larger area in the chaotic do- 
main and it starts showing the effect discussed in section 2, 
i.e., it defines preferential directions in the surface of section, 
given by many manifold segments forming narrow bundles 
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in almost parallel directions. Such bundles are clearly visible 
in both Figs. 6a and 7a. 

Let us return to the three projections of the apocentric 
manifolds on the three planes (0,po), (r,p$), (r,6) shown in 
Figs. 6 a, b, c, respectively. The projection of the apocentric 
manifolds on the plane (r,pe) shows an impressively narrow 
correlation between r and pe- For every branch (U, or UU) 
there is a monotonic (on average) dependence of r on pe. 
This is a numerical verification, in the case of the N-Body 
experiment, of the monotonic relation between r and pg, 
that was discussed in subsection 2.3. 

A direct result of this monotonic (on average) depen- 
dence (of the apocentric distances r and the angular mo- 
mentum pe of the orbits on the invariant manifolds) is that 
causes a near-isomorphy between the projection of apocen- 
tric manifolds in phase space (9,p$) (shown in Fig. 6a) and 
in the configuration space (r, ff) (shown in Fig. 6c) . 

A main effect caused by this near-isomorphy is that the 
bundles of preferential directions of the apocentric mani- 
folds that appear as almost straight lines in Fig. 6a, appear 
as spiral lines in Figs. 6c. In other words the loci of prefer- 
ential directions defined by the apocentric manifolds define a 
geometric pattern on the configuration space that, except for 
some details, has the general form of a grand design spiral. 

It should be stressed that all the points of the apocentric 
manifolds in Figs. 6, 7a,b correspond to local apocentric 
passages of the respective orbits. Thus, one single orbit does 
not describe continually in time the manifolds shown in these 
figures. But if we take a chain of stars that densely populate 
the manifolds, this chain defines a flow along the manifold 
that keeps the apocenters of chaotic orbits correlated for 
quite long times. This correlation was modelled theoretically 
in subsection (2.2) as a soliton-type flow of the successive 
positions of the apocenters of chaotic orbits. 

In Fig. 8 the positions of maxima of the particle distri- 
bution along the N-body spiral arms (Fig. 8a, same as Fig. 
5a) are superposed on the spiral patterns generated by the 
preferential directions of the unstable invariant manifolds of 
the Pli,2 family (plotted for three different values of the Ja- 
cobi constant in Figs. 8b, c, d, which are in a narrow band 
of values larger than the Jacobi constant at the Lagrangian 
points Li,I/2, and close to the peak of the distribution of the 
energies of Fig. 5b). It is immediately seen that the mani- 
fold spirals are very closely aligned to the real spiral arms. 
The images of the apocentric manifolds on the configura- 
tion space, for the three different values of Ej, are similar 
to each other, implying that their superposition enhances 
the overall spiral pattern induced by the manifolds. 

The fact that the projections of apocentric manifolds in 
the configuration space are well-correlated with the maxima 
of the spiral density is remarkable, because the form of the 
manifolds reflects topological properties of the underlying 
chaotic phase space, while the density maxima are related 
to the particles' (chaotic) motion in configuration space. We 
can make a number of comments pointing to the theoretical 
understanding of this phenomenon: 

i) We have compared the distribution of all the parti- 
cles in chaotic orbits of the system along the Jacobi con- 
stant axis given in Voglis et al. 2006 (Fig. 15 therein) with 
the distribution of the particles shown in Fig. 5b. The two 
distributions are quite similar, and their maxima are pro- 
jected at about the same value of Ej, which is very close to 



the band of values of Ej for which the manifolds of Fig. 8 
are calculated. This means that the backbone of the spiral 
arms is related to chaotic orbits with initial conditions close 
to the apocentric manifolds of figure 8. Particles with values 
of the Jacobi integrals much different than those near the 
maximum of the distribution of Fig. 5b partly contribute to 
the density along the spiral arms. However, we have checked 
that most of these particles contribute to the density of the 
axisymmetric background. 

ii) Near the bundles of preferential directions, the man- 
ifolds act, across these directions, as attractors of the orbits, 
since any element of an area-preserving phase flow has to 
contract in the direction across the manifolds in order to 
compensate for the strong stretching along the manifolds. 
In fact, the stretching and twist properties of the manifolds 
are far from being uniform along their length. The attraction 
across the manifold, as well as the organization of phases, 
is stronger where the stretching along the manifold is also 
stronger. This mechanism is more effective along bundles of 
preferential directions, while it is less effective near the turn- 
ing points of the manifold. These phenomena are similar to 
the 'stickiness' phenomena caused by invariant manifolds in 
model Hamiltonian dynamical systems (e.g. Efthymiopoulos 
et al. 1997, Contopoulos et al. 1999). 

iii) Finally, another property of the manifolds is that 
they are recurrent (Contopoulos and Polymilis 1993), i.e., 
the manifold lobes return many times near the points Pl 1 2 
even if they temporarily go to large excursions away from 
Pl-l 2 ■ This mechanism enhances further the density along 
the segments of strong stretching, and also agrees with the 
maxima of the density appearing preferably along bundles of 
the manifold. The recurrence of the manifolds is particularly 
important in the system under study, because the phase- 
space is non-compact and the recurrence of chaotic orbits 
cannot be established by the Poincare recurrence theorem 
which is applicable only to compact systems. It is remark- 
able that the recurrence of the manifolds can create new 
stickiness effects (recurrent stickiness), leading to further ac- 
cumulation of points in the sticky region (Contopoulos and 
Harsoula 2006). 

Figure 9 shows a comparison of invariant apocentric 
manifolds with real spirals at three different time snapshots 
of the evolution of the N-Body experiment, namely t = 18 
(Figs. 9a,b), t = 47 (Figs. 9c,d) and t = 74 (Figs. 9e,f). (The 
manifolds are plotted after eight iterations, when the spread 
of numerical points along the manifold is still small enough 
to allow for a coherent visualization of the manifold). 

The main difference in the system at these three snap- 
shots is that, as t increases, the pattern speed £l p of the 
bar has a tendency to decrease as a result of torques ex- 
erted to the bar by the rest of the system and the angular 
momentum transference outwards. The decrease of Q p is sig- 
nificant, of order 20%, between t = 18 (fi p = 0.8037) and 
t — 47 (Q p — 0.6662), while it is less important (about 6%) 
between t = 47 and t = 74 (Q p = 0.6282). 

As the system slows down, the position of Z/1,2 moves 
outwards and the Jacobi constant -Ej,l 1>2 increases. The un- 
stable manifolds in Figs. 9b,d,f are calculated for the peri- 
odic orbit of the corresponding Pli,2 family at the Jacobi 
constant values Ej = -1.25 x 10 6 , Ej = -1.131 x 10 6 
and Ej = —1.085 x 10 6 , respectively. These values are 
slightly larger than the Jacobi constant values at the re- 
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spective Lagrangian points Li,2, namely Ej — —1.26 x 10 6 , 
Ej = -1.133 x 10 6 and Ej = -1.088 x 10 6 . 

Clearly, as the pattern speed of the bar decreases, the 
manifolds become more and more open, that is they ex- 
tend to larger distances beyond corotation under the same 
number of iterations. Furthermore, Figs. 9b,d,f show that 
the agreement of the apocentric manifold with the spiral 
pattern is good precisely in the regions where the manifold 
segments accumulate in bundles of preferential directions. 
The overall connection of the manifold's spiral geometries is 
generally visible in all three figures, except for a weak spiral 
extension in Fig. 9a corresponding to the outermost points 
of the spiral maxima of Fig. 9b (this is probably a transient 
effect). 

We conclude that as the manifolds become more and 
more open, the bundles of preferential directions become 
broader, but also the real spiral becomes weaker and less 
clearly defined. A detailed explanation of these phenomena 
is necessary, since they all point to the direction of a clear 
connection between the dynamics induced by invariant man- 
ifolds and the dynamics of real spiral arms. In particular, 
given that the locus of apocenters defined by the manifolds 
is the place where stars spend a large part of their radial 
period, we conclude that the projections of the apocentric 
manifolds reveal a narrow dynamical connection between 
manifolds and spiral arms. 

Finally, Fig. 10 shows a calculation of a number 300 of 
individual stellar orbits (plus another 300 symmetric orbits) 
with initial conditions along the unstable direction of the 
manifolds of Fig. 8b. The twelve different panels correspond 
to twelve successive positions of the pericenters (white thick 
dots) and apocenters (black thick dots) of these orbits. The 
first panel (left-top) is when the apocenters start leaving the 
neighborhood of Pl\,li- This figure summarizes the previ- 
ous conclusions from the point of view of the orbital flow. 
Namely, we observe that the phases of these orbits remain 
very well correlated for nine radial periods, and to some 
extent for another three periods. Furthermore, apocenters 
are preferably aligned along the spiral pattern, while peri- 
centers are located mostly in a double ring-like locus sur- 
rounding the bar and the stable Lagrangian points L4, L 5 . 
It is remarkable that inside the corotation radius pericen- 
ters and apocenters cooperate to support a shell of chaotic 
orbits surrounding the bar. Outside the corotation radius, 
for some azimuthal extent from the Lagrangian points L\, 
L2, pericenters and apocenters also cooperate to form an 
enhanced initial arc of the spiral arm. For larger azimuthal 
angles, however, the pericenters form an arc surrounding the 
stable Lagrangian points L4, L 5 , while the apocenters form 
a more open arc following the spiral arms. 

Thus, we conclude that the apocentric manifolds are 
more important than pericentric manifolds, for the forma- 
tion of spiral arms, not only because of the longer time in- 
tervals spent by the orbits near the apocenters than near the 
pericenters, but also because only the apocentric manifold 
is responsible for the extension of the spiral arms to large 
distances from the corotation radius. 



4 CONCLUSIONS 

We have explored the role of the invariant manifolds ema- 
nating from simple unstable periodic orbits (forming fam- 
ilies terminating at L\ ,or L2) in producing a mechanism 
of phase correlations between chaotic orbits beyond corota- 
tion in a barred spiral galaxy. In particular, we define the 
apocentric (or pericentric) one- dimensional manifolds as the 
intersections of the phase space manifolds by the surface of 
section of the apocenters (or the pericenters). A theoreti- 
cal analysis is made, based on a theory of soliton-type mo- 
tions of phases (Voglis 2003) of orbits with initial conditions 
along the apocentric manifolds. A numerical analysis is also 
carried out, referring to the calculation of invariant apocen- 
tric manifolds in the potential of an N-Body simulation of a 
barred spiral galaxy. In this case we explore the connection 
between the manifolds and the overall dynamical features of 
the N-Body system. Our main conclusions are as follows: 

1) While every orbit with initial conditions on the un- 
stable apocentric manifolds of the Pl\,i unstable periodic 
orbits is chaotic, the phases (the positions of the apocenters) 
for a chain of particles with initial conditions populating the 
manifold remain correlated for quite long times. This corre- 
lation can be approximated by a soliton-type motion for the 
phases. This is true not only for small amplitude of the per- 
turbation, but also, to a large extent, for large amplitude of 
the perturbation. 

2) The portraits of phase space, as described by the 
apocentric manifolds, have nearly-isomorphic portraits into 
configuration space. These portraits are described in sub- 
section 2.3 as well as in Fig. 6. In particular, for small val- 
ues of the perturbation, the apocentric manifolds form a 
double-ring structure surrounding the bar and the stable La- 
grangian points Li,Lz- But when the perturbation is strong 
enough, the images of the manifolds in configuration space 
extend to large distances beyond corotation, and they take 
the form of trailing spiral arms. 

3) A comparison of unstable invariant apocentric man- 
ifolds with the density maxima of an N-Body experiment 
along the spiral shows an overall connection of the manifold 
spirals with the N-Body spirals. In particular, the two sets 
are closely aligned in all regions where the manifold seg- 
ments form narrow bundles along preferential directions of 
strong stretching. 

4) As long as the galaxy's pattern slows down, the man- 
ifolds become gradually more open and the bundles of pref- 
erential directions become broader, but still the connection 
between manifold and spiral geometry is preserved. On the 
other hand, the amplitude of the spiral becomes also weaker, 
and its shape less clearly defined. In conclusion, we find a 
clear dynamical connection between the loci of apocenters 
defined by the manifolds of unstable periodic orbits and spi- 
ral arms. In this respect, the apocentric manifold is narrowly 
associated with the positions of maxima of the density along 
the spiral arms. Furthermore, this connection is preserved 
for an important fraction of the simulated galaxies' lifetimes. 

As a final remark, we may revisit the question posed in 
the introduction on whether the above mechanism can help 
clarify the role of the bar in driving the dynamics of spiral 
arms. Our findings suggest that, through the manifolds, the 
bar can indeed initiate a spiral structure. Notice that the 
unstable manifolds are asymmetric in configuration space 
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with respect to the bar axis (fig. 4) and, as shown in section 
2, they can initiate trailing spirals, provided that the bar is 
strong and/or rotates fast enough to create a large effective 
perturbation near corotation. Furthermore, this mechanism 
exploits the relative flexibility of chaotic orbits to exhibit 
large excursions in phase space while remaining correlated 
as regards their phases. On the other hand, a spiral pattern 
initiated by a bar through the above mechanism amplifies 
the effective perturbation even more, a fact that acts addi- 
tively to the whole phenomenon, which is based on chaotic 
orbits. These points, which may prove to be crucial for the 
longevity of the induced spiral pattern, are a subject of cur- 
rent research. 

Acknowledgements: This project was supported in 
part by the Empeirikion Foundation. Stimulating discus- 
sions, as well as a very careful reading of the manuscript 
by Prof. G. Contopoulos are gratefully acknowledged. 



Sellwood, J. A.: 2000, Astrophysics and Space Science 272, 
31. 

Smale, S: 1967, Bull. Am. Math. Soc. 73, 747. 

Sparke, L.S. and Sellwood, J.A.: 1987, Mon. Not. R. Astr. 

Soc. 225, 653. 

Voglis, N., Contopoulos, G., and Efthymiopoulos, C, 1998, 
Phys. Rev. E, 57, 372. 

Voglis, N., Contopoulos, G., and Efthymiopoulos, C, 1999, 
Cel. Mech. Dyn. Astr., 73, 211. 

Voglis, N., Kalapotharakos, C, and Stavropoulos, G., 2002, 
Mon. Not. R. Astr. Soc, 337, 619. 
Voglis. N.: 2003, Mon. Not. R. Astr. Soc. 344, 575. 
Voglis, N., Stavropoulos, G., and Kalapotharakos, C, 2006, 
MNRAS, (in press). 



REFERENCES 

Allen, A.J., Palmer, P.L., and Papaloizou, J, 1990, Mon. 
Not. R. Astr. Soc. 242, 576. 

Athanassoula, E.: 1980, Astron. Astrophys. 88, 184. 
Contopoulos, G.: 1978, Astron. Astrophys. 64, 323. 
Contopoulos, G.: 1985, Comments on Astrophysics 11, 1. 
Contopoulos, G. and Grosb0l, P.: 1986, Astron. Astrophys. 
155, 11. 

Contopoulos, G. and Grosb0l, P.: 1988, Astron. Astrophys. 
197, 83. 

Contopoulos, G. and Polymilis, C: 1993, Phys. Rev. E47, 
1546. 

Contopoulos, G and Harsoula, M. : 2006 (in preparation). 
Contopoulos, G., Harsoula, M., Voglis, N., and Dvorak, R.: 
1999, J. Phys. A 32, 5213. 

Contopoulos, G.: 2004, "Order and Chaos in Dynamical 
Astronomy" , Springer- Verlag. 

Goldreich, P. and Tremaine, S.: 1978, Astrophys. J. 222, 
850. 

Efthimiopoulos, Ch., Contopoulos, G., Voglis, N. and Dvo- 
rak, R.: 1997, J. Phys. A: Math, and Gen. 30, 8167. 
Grobman, D.M., 1959: Dokl. Akad. Nauk SSSR 128, 880. 
Hartman, P., 1960: Proc. Am. Math. Soc. 11, 610. 
Kaufmann, D.E. and Contopoulos, G.: 1996, Astron. As- 
trophys. 309, 381. 

Lynden-Bell, D. and Barot, J.M.: 2003, in G. Contopoulos 
and N. Voglis (Eds) "Galaxies And Chaos", Lect. Notes. 
Phys. 626, Springer- Verlag, pp30-43. 

Patsis, P.A., Contopoulos, G. and Grosb0l, P.: 1991, As- 
tron. Astrophys. 243, 373. 

Patsis, PA., Hiotelis, N., Contopoulos, G. and Grosb0l, P.: 
1994, Astron. Astrophys. 286, 46. 

Patsis, PA. and Kaufmann, D.E.: 1999, Astron. Astrophys. 
352, 469. 

Peebles, P.J.E., 1969, Astrophys. J 155, 393. 

Poincare, H., 1892, "Les Methodes Nouvelles de la 

Mecanique Celeste", Gauthier-Vilars, Paris. 

Romero Gomez, M., Masdemont, J.J., Athanassoula, E. 

and Garcia Gomez, C: 2006, |astro-ph/0603124] 

Schwarz, M.P.: 1984, Mon. Not. R. Astr. Soc. 209, 93. 

Schwarz, M.P.: 1985, Mon. Not. R. Astr. Soc. 212, 677. 



Invariant manifolds and spiral arms in galaxies 11 




Figure 1. Short period orbits around the Langrangian equilibria Li,2 (unstable) and £4,5 (stable). The curves with arrow correspond 
to the equipotcntial passing from Li,2- The arrows indicate the sense of rotation inside and outside corotation. 
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Figure 2. The Poincare surface of section {12,82) for 6\ = 2kn,k = 0,1,2,... in the Hamiltonian (0 with k* = 3.23414 X 10 _1 , 
o* = -2.90577 X 10" 2 , 6* = 0, c = -7.5 X 10" 3 , A, = 1.897 X 10~ 2 and (a) e = 10~ 3 , (c) e = 10~ 2 . In both cases the value of the 
Jacobi constant is set equal to Ej = —1.897 X 10~ 2 . The stable and unstable manifolds of the periodic orbits Pli,2 (unstable fixed 
points in (a),(c)) are shown in (b) and (d) respectively. 
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Figure 3. Comparison of the theoretical soliton-type solution given by Eq.([l3}, (solid line) with the corresponding numerical relation 
for 200 points along the unstable manifold U emanating from Pjj\ in the Hamiltonian (|3]l with parameters as in Fig. 2d. The different 
panels correspond to different time snapshots, namely a)t = 0, b) t = 120, c) t = 240. 
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Figure 4. Schematic representation of the isomorphy between a) the phase portrait in the Poincare section (I2, O2), and b) configuration 
space. The regions A,B,C,D in (a) correspond to A', B', C, D' in (b). The arrows indicate the directions of the stable and unstable 
manifolds emanating from Pli,2 in phase space and configuration space. The manifolds emanating from Pli 2 are shown in phase space 
(c,e), and configuration space (d,f) for a short time (four iterations of an initial segment ds = 10~ 4 ), (c,d) and longer time (12 iterations), 
(e,f). 
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Figure 5. a) Projection of the particles of the N-Body simulation described in section 3 on the (y, z) plane, at the time snapshot t = 47. 
In this projection the galaxy rotates counterclockwise. The thick dots are the maxima of the surface density calculated with a division of 
the system in concentric rings (see text for details), b) The distribution of the Jacobi constant values (in units of the N-body code) for 
the particles located in small bins ArAS centered at the thick dots of (a). The dashed vertical line indicates the Jacobi constant value 
for L\. 
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Figure 6. Projections of the apocentric manifolds: (a) on the surface of section plane (9,pg), (b) on the plane (r,p$) and (c) on the 
configuration space (r, 9) or (z,y). In (a) and (b) the unstable fixed point Pl2 and the corresponding initial directions of the manifolds 
U and UU are also shown. In (c) the projection of the symmetric manifolds, emanating from the fixed point Pli, is shown together with 
the manifolds from Pli. The monotonic (on average) dependence of r on pg shown in (b) produces a ncar-isomorphy between (a) and 
(c). 
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Figure 7. a) Unstable invariant manifolds of the unstable periodic orbit P^2 in the Poincare section (8,pg), f = 0, r < of the potential 
of the N-Body simulation, at the time snapshot t = 47. The initial directions are marked U and UU. The initial segments along U and 
UU, of length ds = 1Q~ 4 , are populated with 10000 points each. The figure shows the image of these points after six iterations of the 
Poincare map. The Jacobi constant is Ej = —1.125 X 10 6 . b) The manifolds of (a) projected in configuration space, together with the 
symmetric manifolds with respect to the center, c) The evolution r vs. r for two orbits which give the points of intersection P and P' in 
the Poincare section, d) The time evolution of r(t) for the same pair of orbits. 
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Figure 8. a) Same as in figure 5a. The form of the spiral (thick dots) is compared with the forms of the unstable manifolds emanating from 
the unstable periodic orbits Pli,2 f° r three different values of the Jacobi constant, namely b) Ej = —1.131 X 10 6 , c) Ej = —1.128 X 10 6 , 
and d) Ej = -1.125 x 10 6 . 
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Figure 9. The projection of the N-Body system on the (y, z) plane at three different snapshots: a) t = 18, c) t = 47, and e) t = 74. The 
comparison of the spiral pattern with the unstable invariant manifolds emanating from Pli,2 is made in (b), (d) and (f) respectively. 
The corresponding Jacobi constant values at which the manifolds arc calculated is Ej = -1.25 X 10 6 for (b), Ej = -1.131 X 10 6 for (d) 
and Ej = -1.085 X 10 6 for (f). The 8-th manifolds' iteration is shown. 
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Figure 10. Twelve successive Poincare consequents of the apocenters (thick black dots) and the pcricenters (thick white dots) of 600 
orbits with initial conditions on the unstable manifolds U emanating from the fixed points Pli,2 (300 orbits for each fixed point). The 
first panel corresponds to the Poincare consequent at the time when the apocenters start leaving the neighborhood of the fixed points. 
The phases of the apocenters of these orbits remain well correlated for at least 9 successive passages. The phases of the pericenters are 
well correlated for much longer. The pericenters tend to form a double-ring surrounding the bar and the stable Lagrangian points L4, 
L5 and cooperate with apocenters to support a chaotic shell surrounding the bar, but also an azimuthal arc of the spiral arms near L\ , 
L2, outside corotation. Spiral arms at larger azimuthal angles and larger radii can be supported only by apocenters. 



